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Abstract 

The classically perfect action of QCD requires no tuning to get the pion 
massless in the broken phase: the critical bare mass mj is zero. Neither the 
vector nor the flavour non-singlet axial vector currents need renormalization. 
Further, there is no mixing between four-fermion operators in different chiral 
representations. The order parameter of chiral symmetry requires, however 
a subtraction which is given here explicitly. These results are based on the 
fact that the fixed point action satisfies the Ginsparg- Wilson remnant chiral 
symmetry condition. On chiral symmetry related questions any other local 
solution of this condition will produce similar results. 



1 Work supported in part by Schweizerischer Nationalfonds. 



1 Introduction 



Lattice regularized, local QCD actions break chiral symmetry under general con- 
ditions Wilson fermions Q have a dimension five symmetry breaking operator 
whose effect on the physical predictions goes to zero in the continuum limit. It 
leaves, however its trace behind in the form of additive quark mass renormalization, 
axial current renormalization and mixing between operators in nominally different 
chiral representations. There is a significant recent progress in calculating these 
renormalizations in a theoretically controlled non-perturbative way Neverthe- 
less, the situation is not really pleasing theoretically and the technical difficulties 
are also significant. 

Staggered fermions Q keep a part of chiral symmetry intact which offers an- 
other possibility to study problems where chiral symmetry is essential ||. On the 
other hand, staggered fermions do not solve the doubling problem and constructing 
operators with correct quantum numbers is far from trivial. 

A rather different method to overcome the problems with chiral symmetry is 
the domain wall fermion || or the overlap formalism Q . Recent results on kaon 
matrixelements are nice and promising 0. Similarly to the staggered fermions, the 
domain wall fermions were motivated by a single issue: to solve the problems of 
chiral symmetry in the fermion sector. 

This paper is part of a project to construct a lattice formulation for QCD which 
performs well in every respect both in the gauge and in the fermion sectors, includ- 
ing classical solutions, topology, cut-off effects and chiral symmetry. The fixed point 
(FP) action (m q = 0) and the actions on the trajectory along the mass direction 
(m q 0), which are local and determined by saddle point equations, are classically 
perfect [|[ [k| |l^|. Among others, the FP action is perfect concerning the clas- 
sical solutions leading to scale invariant instantons |)| [To[ [l3j [ll| and fermionic 
chiral zero modes satisfying the index theorem fl6|| . Even more, as we discuss here, 
these actions are quantum perfect what concerns chiral symmetry |jl7| . The problem 
of constructing the FP action and the mass trajectory by solving the corresponding 
classical saddle-point equations requires skill and patient. The real difficulty is, 
however to find a parametrization for these actions which is sufficiently precise and, 
at the same time, does not make the simulations too expensive. Examples in d = 2 
, [l8| that with a limited number of couplings a parametrization can be 
achieved which performs excellently in simulations and other numerical checks. In 
d=4, although the preliminary results look promising, the parametrizations studied 
until now are admittedly rather primitive [ fl0[ ^9[ . 

As the discussion above shows, there exist several possibilities to cope with the 
problems of chiral symmery. It might turn out on the long run that none of the new 
theoretical ideas can compete with the direct approach of using Wilson fermions 
and calculating renormalization factors, mixing coefficients in different processes 
and driving the system close to the continuum limit to kill the cut-off effects and 
to avoid other problems ^0|. The author hopes that the solution will be more 
appealing. 

It is somewhat surprising that the FP action respects chiral symmetry in its 
predictions, since - complying with the Nielsen-Ninomiya theorem 0] -, it breaks 
chiral symmetry explicitly. The basic observation is that this breaking is realized 
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in a very special way. The FP action satisfies the equation |T^, [T(| |^| 

\{Kn>(UWW = Rnn>{U), (1) 

or equivalently 

i{/»n„',7 5 } = (^7 5 ^W. (2) 

Here /i(C/) — 1 is the FP quark propagator over the background gauge field U and 
the matrix R nn ' is trivial in Dirac space. Both h(U) nn ' and R{U) nn i are local. The 
precise form of R depends on the block transformation whose FP we are considering. 

It has been observed a long time ago that eq. ( 3) w ith a local R is the mildest 
way a local lattice action can break chiral symmetry pj]|. Along the mass trajectory 
(m q 0) eq. (||) is modified as 

-^{Ln', 7 5 } = (h 5 Rh)nn> + m,I& n , , (3) 

where we denoted the action along the mass trajectory by h, h — > h as m q — » 0. 
The operator E 5 is a local pseudoscalar density which goes to 7 5 <5 nrl ' in the formal 
continuum limit. The results of this paper are based on eq. (|^), hence any local 
solution of this equation will lead to the same results. It is interesting in this context 
that in a recent paper p2|| another solution (growing out of the overlap formalism) 
of eq. ([!]), which seems to be unrelated to renormalization group considerations, was 
presented. This is an interesting development, if it can be shown that the Dirac 
operator in p^] is local over any gauge configuration. 

Since the locality of h, R and the currents is basically important in deriving the 
following results, let us discuss briefly what is meant by this notion. An operator 
(like action density, topological charge density, current, etc.) on the lattice is local 
if it has an extension of 0(a): the coupling between the fields in the operator at 
a distance n decays exponentially as exp(— jn), with 7 = 0(1), or it is identically 
zero beyond a certain O(a) range. In the continuum limit (a — > 0) the extension 
measured in physical units goes to zero. A quantum field theory with non-local 
interactions looses renormalizability and universality in general - this would be 
certainly too high a price for an 'improvement'. 

The remark has been made already in |21j that the soft-pion theorems are ex- 
pected to remain valid if eq. (|2|) is satisfied. The intuitive reason is that the two 
h factors on the r.h.s. of eq. (|2|) will cancel the two propagators which connect 
this term to other operators in the matrixelement producing only a contact term in 
Ward identities. On the other hand, in current algebra relations, in Ward identities 
related to mixing and elsewhere, contact terms are relevant, so it is necessary to 
study in detail what eq. (^|) implies on chiral symmetry. Many of our considera- 
tions rely on the seminal paper of Bochicchio et al. [ p3| where the chiral symmetry 
properties of Wilson fermions are discussed. 



2 We denote the commutator and anticommutator by [,] and { , }, respectively. 
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2 The vector and axial currents 



The classically perfect fermion action has the form 

Sf P ($, V, U) = YJ Tpm>h m >n>(U)lp n > . (4) 
m' ,n' 

The Dirac, flavour and colour indices are suppressed in eq. (Q) . The Dirac operator 
h is assumed to be flavour independent. The flavour generators are denoted by 
r a ,a=l,...,Nj-l, tr{T a T b ) = \5 ab . 

As usual, the axial currents are constructed from the chiral symmetric part of 
the action using /isym = ^[h,j 5 ]j 5 . There is a certain freedom in constructing the 
vector currents. Using the full action, conserved vector currents are obtained. On 
the other hand, the asymmetry between the vector and axial case generates strange 
terms in current algebra relations and it requires extra work to show that they go 
away. We decided to follow and construct both currents with /isym- We get 

V„V£(n) = $ n r a (hsYMip) ~ (<Msym) r a V„ , (5) 



V^(n) = -$ n T a (7%™^) ~ (^sym7 5 ) n r a ipn ■ (6) 

None of these currents are conserved, at least not without further considerations, if 
the equations of motion hip = is used. 

In order to find the currents themselves we introduce the flavour gauge matrices 
W^in) = 1 + iw°(n)r° + . . . and extend the product of colour U matrices in eq. (|[) 
along the paths between the fermion offsets m',n' by the product of W matrices 
along the same paths. Then the vector current is defined by 



^"W = -i 53 VW Swa ( n ) (WM)m'»' (U, W)l/jr, 



(7) 



This current satisfies eq. (0) Jm) . Every path in f /isym ) which goes through 

V / m'n' 

the link (n, n + jl) gives a contribution to V£(n) which is equal to the contribution to 
( /isym ) times (s+ — s_)r a , where s+(s_) is the number this path runs through 

V / m'n' _ 

(n, n + (i) in positive (negative) direction. Eq. (j7|) can also be written as 

V£(n) = n; U)i> n . , (8) 

m'n' 

where T n is trivial in flavour space. Similar considerations lead to the axial current 
A a a (n) = ^m'T a Tl(m', n'; n; U)^ n , , (9) 

m'n' 

where 

it = r P7 5 = -7% . (io) 

3 We use the standard sign convention for the currents as opposed to that in [l6| . We use 
the notation V M /(n) = f(n) — f(n — fl) . 
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3 The basic Ward identity 



We shall repeatedly use the following Ward identity which can be obtained by 
changing integration variables in the QCD path integral (see, for example E3j) 



i(— 0{ Xl ,..., Xj )) = {0(x 1 ,...,x j )V ll Al(n)) 



S_ 

(0(x 1 ,...,x j )X a (n))- (11) 
(0(x 1 ,...,x j )2m q P a (n)), 

where O is some product of local operators, SO is its change under a local chiral 
transformation with the infinitesimal parameter e°, P a (n) is a pseudoscalar density 

P a (n) = \$nT a (S 5 V)„ + (V^ 5 )„t>„] , (12) 
while X a comes from the remnant chiral symmetry condition eq. (g) 

X a (n) = ^{h^Rh^n + (fh^Rh^T^n . (13) 

In deriving eq. ( [TT|) we used eq. (^), hsYM = h — ^{h,j 5 }^ 5 and eq. (||). For later 
use we introduce the flavour singlet scalar density S(n) 

S(n) = i$„(E^) n + $£) n Vfc] (14) 



with 



S = ^{S 5 ,7 5 }. (15) 



The limit m q — >• and the order parameter of 
chiral symmetry 



The following considerations are valid in the broken phase in general, not only in 
the continuum limit. We shall consider the Ward identity eq. ( p"l| ) with O = P b (x) 
and sum over n. After some algebra we obtain 

E^ pb(x)) = ~ 5ah ^ (5(x)) ' (16) 



-Y l (P b (^X a (n)) = S ab {tT DC -{E 5 (U), 1 5 R(U)U), (17) 

n 

where the trace tr DC in eq. (|l7|) is over Dirac and colour space. We used the 
notations introduced in eqs. ( |12[]15| ). The last term in eq. ( pd| ) is dominated by the 
pion state for small m q : 

|2 

~2m q {P b (x)J2P a (n)) =2m q 6 ab 1 -^, (m q -> 0) , (18) 

n n 

where we denoted the pion-vacuum matrixclcmcnt of the pseudoscalar density P 
by a. We shall see later (eq. (|l])) that \a\ 2 — /^(m^/m 9 ) 2 /4 in the continuum 
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limit, but it is sufficient at this moment to observe that a is finite as m q — > 0. 
Eqs. give 

1 <rft 

-{S{x)+Nft^ c -{Y?^R} xx ) = 2N f ^L\a\\ (m g - 0) . (19) 

We shall show now that the combination on the l.h.s. is one of the possible order 
parameters of chiral symmetry: it is zero in perturbation theory, or more generally 
in the limit limy—^ lim m(j ^o ( m this order). If this order parameter picks up a 
non-zero expectation value in the limit \im mq ^o limy_>oo due to non-perturbative 
effects then eq. ( |l9| ) implies 

ml ~ m q , (m, -» 0) , (20) 
i.e. the critical quark mass is zero. There is no tuning. 

Let us first construct the simplest order parameter. Consider 

(^x^x)sub = (^xV'x + 4N f tr c R{U) xx ) (21) 
in the limit m q — ► in a finite volume. Integrating out the fermions we get 

(^x)sub = <tr DFC (-^ + R xx )) , (22) 

where the trace is over Dirac, flavour and colour space. The r.h.s. of eq. ( p2[ ) is zero 
as can be seen easily by taking the trace of eq. ([[]) . For specific block transformations 
('blocking out of continuum' plfl) R xx — 1/k, where n is a parameter entering the 
blocking procedure. In this case the order parameter is obtained from (ipxipx) by 
simply adding the constant 4NfN c /K. 

Eq. (|j) allows to construct other order parameters also. Multiplying eq. ([!]) by 
S 5 7 5 from the left and adding to it the product from the right, after taking the 
trace the expectation value on the l.h.s. of eq. ( |l9| ) is obtained as another definition 
of the order parameter. This is what we wanted to show. 



5 Current renormalization 

Current algebra relations which are non-linear in the currents, are very convenient 
for the study the renormalization of the currents. Consider the Ward identity, 
studied also by Bochicchio et al. p3[ in this context, in the continuum 

d*(A*(x)A b u (y)V p c (z)) = 2m q (P a (x)AUy)VC(z)) + 

if abd S(x-y)(V v d (y)V;(z)}+ (23) 
if acd 8{x-z){A b v {y)A%z)). 

Integrating over x, we obtain 

= J dx2m q (P a (x)Al( v )V^z)) + if M (V u d (y)V p c (z)) 

+ ir d (A b u {y)A d p {z)). (24) 

We shall consider eq. ( |2~4| ) for \y — z\ being a physical distance. With this condition 
we assure that the communication between the points y and z goes through physical 



G 



intermediate states. We work out the Ward identity corresponding to eq. ( f24| ) on 
the lattice. Consider the general Ward identity in eq. (TO) with 



O^Al(y)V p c (z), (25) 
where \y — z\ is much larger than the lattice unit a. We get ^ 



J2i^Al(y)=if abd V v d (y). (27) 



Using eqs. d25|j2qj27D in eq. ([Ill), the continuum Ward identity eq. fl24| ) is repro- 
duced if 

{J2x a (n)A b MV p c (z)) = 0. (28) 

n 

Eqs. (pPjnf) give after integrating out the fermions 

- 2Tr (VrVXiyzh^lt (y) + Tl(y)j 5 RT p {z))%- x ) , (29) 

where Tr is a trace over all the indices, including space. Since R, T p and are 
local, eq. ( p9| ) is zero for \y — z\ much larger than the lattice unit. We get, therefore 

= £2m 9 (P>)^(2/)V p c (z)) + if abd (V d (y)V p e (z)) + ir d {A b v {y)A d p {z)) , 

(30) 

where P a (n) is defined in eq. ( p"2| ) . 

Consider finally the Ward identity eq. ( |Tl| ) with O(x) — P b (x) as in Section 4 
but do not sum over n. Assume that x — n is much larger than the lattice unit. 
Using similar steps as above we obtain 

(P 6 (z)V^») - (P b {x)2m q P a {n)) . (31) 

From eq. ([}!) we conclude that V^^l" and 2m q P a (n) have the same renormalization 



factor. Eq. ( |30| ) gives then 

Z A = Z v = 1 . (32) 
We remark that in the considerations above the m q — > limit was not necessary. 

6 Mixing 



In the study of weak, non-leptonic matrixelements it is essentially important to 
construct operators in definite chiral representations. In the case of Wilson fermions 
the nominal (tree level) chiral assignment of operators is invalidated by quantum 

4 The simple relations in eqs. (MJ27|) would receive additional terms had we used the conserved 
form of the vector current. Although they can be shown not to influence the final conclusions, 
their presence would complicate our considerations. 
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corrections due to the chiral symmetry breaking terms in the action. Operators 
with definite chiral properties are linear combinations of the operators with nominal 
assignment (mixing). For present Monte Carlo calculations, as the case of the Bk 
parameter shows, the mixing coefficients should be calculated non-perturbatively, 
which is a highly non-trivial problem. 

Using the example of local k-fermion operators (k=4,6, ... ) we are going to 
demonstrate now that k-fermion operators in a definite chiral representation (assign- 
ment on the tree level) do not mix with other k-fermion operators from a different 
chiral representation. In particular, for the As = 2, 4-fermion operator, whose 
matrixelement between K° and K° defines the Bk parameter, the nominal chiral 
assignment is equal to the full quantum chiral assignment - there is no mixing. 

Consider the Ward identity eq. ( |TT|) with O — > 0(x)B(yi, . . . ,Dj), where 0(x) 
is some k-fermion operator in a definite chiral representation, B{jj\ 1 . . . , yj) is some 
product of local operators, and \y\ — x\, .. . , \yj — x\ are much larger than the lat- 
tice unit. Choosing the point n in eq. (|ll|) equal to x, or, if the operator O has 
an extension of 0(a), summing over n in the neighbourhood of x, then on the 
l.h.s. the chiral variation of the k-fermion operator 0(x) enters. This is a linear 
combination of k-fermion operators in the same representation. On the r.h.s., the 
(X a (n)0(x)B(yi, . . . ,yj)) matrixelement might also produce k-fermion terms local 
in x — n in different representations. This is the source of mixing when using Wilson 
fermions p3| . In our case 

(X a {n)0{x)B{ yil ... lV3 ))= (33) 
((^nT a (hj 5 RhiP) n + {fh^Rh) n T a ^ 0{x)B( yi , . . .,%•)) . 

If one of the fermions in X a are paired with one of the fermions in B, the propagator 
is cancelled by h in X a and the result will be zero since n ~ x is far (in lattice 
units) from y\,...,yj and R and h are local. Hence, both fermions in X a {n) should 
be paired with the fermions in the k-fermion operator 0{x). The result is a (k- 
2)-fermion operator local in the point x. The dangerous term, therefore can not 
give a k-fermion contribution to the l.h.s. of eq. ([fl]). In the example of the Bk 
parameter, the matrixelement in eq. ( |33| ) is zero, since no 2-fermion operator exists 
with As = 2. 
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